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Our aim is to explore the relation between gas, atomic and molecular, and dust in spiral galaxies. Gas surface 
densities are from atomic hydrogen and CO line emission maps. To estimate the dust content, we use the disk 
opacity as inferred from the number of distant galaxies identified in twelve HST/WFPC2 fields of ten nearby 
spiral galaxies. The observed number of distant galaxies is calibrated for source confusion and crowding 
with artificial galaxy counts and here we verify our results with sub-mm surface brightnesses from archival 
Herschel-SPIRE data. We find that the opacity of the spiral disk does not correlate well with the surface 
density of atomic (Hi) or molecular hydrogen (H2) alone implying that dust is not only associated with 
the molecular clouds but also the diffuse atomic disk in these galaxies. Our result is a typical dust-to-gas 
ratio of 0.04, with some evidence that this ratio declines with galactocentric radius, consistent with recent 
Herschel results. We discuss the possible causes of this high dust-to-gas ratio; an over-estimate of the dust 
surface-density, an under-estimate of the molecular hydrogen density from CO maps or a combination of 
both. We note that while our value of the mean dust-to-gas ratio is high, it is consistent with the metallicity 
at the measured radii if one assumes the Pilyugin & Thuan calibration of gas metallicity. 
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1 Introduction 

The radio 21-cm emission of atomic hydrogen (Hi) ob- 
served in the disks of spiral galaxies is a powerful tracer 
of the presence and dynamics of the interstellar medium 
(ISM) , extending to well outside the typical scale of the 
stellar disk. Its origin is likely a mix of "primordial" 
(Fall & Efstathioul 1980), or recently accreted mate- 



rial (Sancisi et al. 2008), recycled matter (ejecta rain- 
ing back onto the disk; e.g., Oosterloo et al. [2007), and 
skins of photo-dissociated material surrounding molec- 
ular clouds ( Allen et al. 2004 1 . The other components 



of the ISM, ionised and molecular hydrogen, metals 
and dust, are all more difficult to trace, because their 
emission strengths depend on the local degree of excita- 
tion which in turn is affected by particle densities and 
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temperatures, photon densities, and stellar and AGN 
illumination. 

Molecular hydrogen is usually traced with CO(J=l- 
or 2-1) line emission, and from it we have derived 
our knowledge of the molecular clouds in nearby spi- 



rals (e.g, Leroy et al. 2008 Rosolowsky 2005). How- 



ever, it remains an open question how sensitive the CO 
brightness is to the local volume density and tempera- 
ture of the ISM, and what is the accuracy with which 
observations of CO surface brightness can be converted 
into H2 column densities and ultimately into molecular 
cloud masses. This conversion is also likely to depend 



on metallicity and hence galactocentric radius (Foyle 



et~al~j [2T7T21 [Israeli [19971 ILeroy et al [ [20071 [2011] [Mad 



den et al. 1997 Pohlen et al. 2010 ) 



Nevertheless, a successful and extensive description 
of the atomic and molecular ISM in spirals and their 
relation to the star-formation rate is currently being de- 
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veloped, using a multi- wavelength approach to estimate 
the star-formation rate, and high-resolution Hi and 
CO observations to characterize the ISM in individ- 



ual galaxies (Bendo et al. 2010b| 


Calzetti et al. 2005 


Foyle et al.||2012| Kennicutt et al. 


2007 Thilkeretal. 


2007 


), in detail in small samples of galaxies ( 


Bigiel 


et al. 


2008 


|Boissier et al.| |2007| |Cortese et al. 2006 


Leroy et al[|2008||Schruba et al. 2011 


), or in a general- 



the dust together with the H I from the disks of spirals 
in a cluster environment. 

In the comparison between the Herschel cold grain 
emission, and H I and CO observations, the mass-opacity 
coefficient of dust grains appears to be too low in M33 
(the inner disk, Braine et al.| 2010[ ), and M99 and 



ized way over a population of galaxies (e.g., Bell et al 



20031 IBuat eTaL] [2002} jCatinella etsl] |2010[ |Fabello| 
[20TT1 |Kannappah] [2004) |Kennicutt| [1998) [West] 
2010). Star- formation occurs when the combined 



et al. 



et al. 



M100 ( |Eales et al.[ |2010j). This is either because (1) 
its value is not well understood, (2) the conversion fac- 
tor between CO and molecular hydrogen, Xco> is dif- 
ferent in M99 and M100, or (3) the emissivity is 
different at sub-mm wavelengths. |Roman-Duval et al. 



ISM exceeds a threshold surface density (although the 
exact threshold is still debated, see e.g., |Bigiel et al. 



2008 Pflamm-Altenburg fc Kroupa 2008). The ratio 



between molecular and neutral ISM is set by the hydro- 



static pressure (jBigiel et al.| |2008[ |Leroy et al.| |2008 



(2010) compare CO, Hi and dust in the Large Mag- 
ellanic Cloud (LMC), and argue that the cause of the 
discrepancy cannot be a different emissivity, nor a dif- 
ferent gas-to-dust ratio, but that CO clouds have H 2 
envelopes, hence Xco changes with different density 



Obreschkow et al.l 120091). Also, observational models environment s (an explanation also favored by |Wolfire 



et al. 2010). Other recent results seem to back varia- 



of the role of photo-dissociation in the balance between 

atomic and molecular hydrogen have made steady progress tlons ln X co5 |Leroy et al.| (|2011| find a link between 



(Allen et al. 


2004 




1997; 


Heiner et al. 


2008a 2009 


2010 


2008b 


Smith et al. 


2000 


)• 



As an alternative to CO, one could use interstellar 
dust as a tracer of the molecular component in spiral 
galaxies, since it is linked mechanically to the molecular 



phase (Allen et al. 1986 Weingartner & Draine 2001 ), 



by mutual shielding from photo-dissociation, and the 
formation of molecular hydrogen on the surface of dust- 
grains (e.g., Cazaux & Tielens 2004). Interstellar dust 



can be traced by its emission or its extinction of starlight 

Surface densities of dust in spirals have been ob- 
tained from spectral energy distribution models of multi- 
wavelength data (e.g., ffioseliretaL 2010 Draine et al 



2007 Popescu et all 12000 Popescu & Tuffs 2002), m 



from simple (modified) blackbody fits of far-infrared 



and sub-mm data ( |Bendo et al.| |2008| |2010b| |Gordon 



et al. 




2008 




2010 


1 or FUV/FIR ratios 


Boissier et al. 


2004 


2007 


2005 




Munoz-Mateos et al. 




20111. The aim 



is to estimate the typical temperature, mass, compo- 
sition and emissivity of the dust, and the implied gas- 
to-dust ratio ( Boissier et al. 120041 iBoselli et al. 



Foyle et al. 2012 



2010 



2011: Munoz-Mateos ct al. 



Galametz et al.[ |2012| |Galliano et al 



Roman-Duval et al. 2010 



2009 a|bl|Pohlen et al.| |20l0 



Smith et al. 2010 ?) 



The most recent Herschel results include a resolved 



temperature gradient in the disks of spirals (Bendo 



et al. 



2012 



2010b Engelbracht et al. 2010 Foyle et al 



Pohlenetal. 2010 Smith et al. 2010 ), linked to 



increased illumination of the grains, notably in the spi- 



ral arms (Bendo et al. 2010b) and bulge (Engelbracht 



et al. 2010). With sufficient spatial sampling, one can 



extract the ISM power spectrum but this is only pos- 



sible with Herschel for local group galaxies (Combes 
et al.[ 2012 ). Based on Herschel data of the Virgo clus- 



ter, (Smith et al., 2010), Cortese et al. (2010) and ? 



show the spatial coincidence and efficiency of stripping 



Xco and metallicity based on SED models of a few 
local group galaxies and the HERACLES CO survey. 
A solid result from the first Herschel observations is 
that the gas-to-dust ratio increases with galactocen- 



tric radius (Pohlen et al. 



do ( |Bendo et al.| |2010a| 



2010 Smith et al. 2010 ) as 



Munoz-Mateos et al. 2009a 



based on Spitzer data alone). ? find a much lower than 
Galactic CO-to-H2 conversion factor based on the re- 
lation between metallicity and gas-to-dust ratio radial 
profiles of several Virgo cluster spirals. 

Even with the excellent wavelength coverage of Her- 
schel, the SED fit results remain degenerate between 
dust mass, temperature and emissivity (see the reviews 
Calzettil 120011 iDrainel 120031). It is still especially 



difficult to distinguish between a mass of very cold 
(poorly illuminated) dust from dust with much differ- 
ent emissivity characteristics (the emissivity efficiency 
depends on wavelength as A - ^ in the sub-mm regime 
with f3 ^ 2, which may be typical for very large grains). 

While large masses of extremely cold dust can be 
ruled out with increasing confidence, the level of illu- 
mination of the grains by the interstellar radiation field 
remains a fully free parameter in the SED models. The 
main uncertainty is complex relative geometry between 
the dusty filamentary structures and the illuminating 
stars. Both the grain emissivity and dust/star geome- 
try can be expected to change significantly throughout 
the disk, i.e., with galactocentric radius or in a spiral 
arm. 

Alternatively to models of dust emission, one can 
use the absorption of stellar light to trace dust densi- 
ties. The advantages are higher spatial resolution of op- 
tical wavelengths and an independence of dust temper- 
ature. However, one needs a known background source 
of stellar light to measure the transparency of a spiral 
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diskQ Two observational techniques have been devel- 
oped to measure the opacity of spirals and consequently 
their dust content. The first one uses occulting galaxy 



2 Galaxy Sample and Data 



pairs 


( Andredakis & van der Kruit[|1992||Berlind et al. 


1997 


Domingue et al. 


1999 


2000; Elmegreen et al. 


2001 


Holwerda et al. 


2007b, 


2009] |Keel et al. 


2012, 


submitted Keel & White| |2001a|b |White et al. 


2000 



Holwerda et al. submitted.), of which an increasing num- 
ber are now known thanks to the Sloan Digital Sky Sur- 



vey and the GalaxyZOO citizen science project (Lintott 



et al. 2008 1 



The second method uses the number of distant galax- 
ies seen through the disk of a nearby face-on spiral, 
preferably in Hubble Space Telescope (HST) images. 
The latter technique is the focus of our "Opacity of 



Spiral Galaxies" series of papers (Gonzalez et al. 1998 


2003 


Holwerda et al. 


2007a 


2005a|b|c|d|e|) | The ben- 



efit of using distant galaxies as the background light 
source is their ubiquity in HST images of nearby galax- 
ies. Now that uniform H I maps are available from the 



THINGS projec t (The Hi Nearby Galaxy Survey, [Wal- 
ter et all] |2008[ ), as well as public Herschel data from 



the KINGFISH (Key Insights on Nearby Galaxies: a 



Far-Infrared Survey with Herschel, 


Dale et al.[ 


Galametz et al. 


2012 


Kennicutt et al. 




2011 ? 



'), and 

CO(J=2-l) mapsTrom the HERAC LES~survey (The 
HERA CO Line Extragalactic Survey, |Leroy et al.||"2009| ) 
for a sub-sample of the galaxies analysed in our "Opac- 
ity of Spiral Galaxies" project, we are taking the op- 
portunity to compare our disk opacities to Hi and H2 
surface densities to see how they relate. 

Our method of determining dust surface densities 
is certainly not without its own uncertainties (notably 
cosmic variance, see but these are not the ones of 
sub-mm emission suffers from (grain emissivity, level of 
stellar illumination, variance within the disk or these). 
Hence, our motivation for our comparison between the 
disk opacity and the other tracers of the cold ISM is 
to serve as an independent check to the new Herschel 
results. 

In section 2, we discuss the origin of our sample and 
data. Section 3 explains how we derive a disk opacity 
from the number of distant galaxies. In section 4, we 
discuss the distant galaxy number as a function of H 1 
column density and in section 5, we compare the Hi 
and H2 column densities, dust extinction, averaged over 
whole WFPC2 fields, and per H 1 contour, respectively. 
Sections 6 and 7 contain our discussion and conclusions. 



Our present sample is the overlap between the [Holw- 



erda et al.l d2005bb, the THINGS (|Walter et al.[p2008), 



and the HERACLES (Leroy et al.fl2009) projects. The 



common 10 disk galaxies are listed in Tabic [T] We use 
the public THINGS data and early science release data 
from HERACLES. 

Figure[TJshows theHST/WFPC2 "footprints" over- 
laid on the VLA HI maps. In the case of NGC 3621 and 
NGC 5194, there are two HST/WFPC2 fields available 
for each galaxy. 

2.1 VLA 21-cm Line Observations 

For this study we use the THING S (The Hi Nearby 
Galaxy Survey, Walter et al. 2008[ ) robustly-weighted 
(RO) integrated total H 1 intensity maps (available from 
|http : //www . mpia- hd . mpg . de/THINGS/ ) . The maps were 
obtained with the VLA, and converted to Hi surface 



density using the prescription from Walter et al. ( 2008 1, 
equations 1 and 5, and Table 3. Although the naturally- 
weighted maps are markedly more sensitive to the largest 
scale H 1 distribution, the robust maps have the highest 
angular resolution. 

The robust maps are better suited for a direct com- 
parison with the number of background galaxies, as we 
are interested in the H 1 column density at the position 
of each background galaxy and hence at scales smaller 
than the FOV of the HST/WFPC2 FOV (3 CCDs of 
l!3x 1.'3). Additionally, we use the WFPC2 footprint 
as an aperture on the Hi maps (Figure [T|[^] The Hi 
column densities averaged over the WFPC2 footprints 
(an angular scale of 2'3) on the sample galaxies, and 
expressed in units of M Q /pc 2 are listed in Table [l] 
These mean column densities include a correction fac- 
tor (1.36) for Helium contribution to the atomic gas 
phase. 



We used the term "opacity" throughout our project and its 
publications for historical reasons. 

2 Other authors have used distant galaxy cou nts or colours 
to estimate extinction in the Magellanic Clouds (IDutra et al. 
200l1|Gurwell fc Hodgel|1990 l |Hod"gel|1974||Hodge fc Snow||1975 
MacGillivrayl |1975| |Shapley[ |1951| |Wesselink| |1961| ) and other 
galaxies I jCuillandre et al.| |2001| |Zaritsky| |1994| . 



3 In this case it does not matter whether the maps are robustly 
weighted or naturally weighted. 
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2.2 CO(J = 2^1) Line Observations 

The HERACLE S project (The HERA CO Line Extra- 
galactic Survey, |Leroy et al.| |2009| |Walter et ah] |2009[ ) 
is a project on the IRAM 30m telescope to map the 
molecular gas over the entire optical disks (R25) of 40 
nearby galaxies via the CO(J=2-l) emission line. The 
HERA instrument has comparable spatial (11") and 
velocity (2.6 km/s) resolutions to the THINGS sur- 
vey, and good sensitivity (3a ~ 3M /pc 2 ) as well. 
The HERACLES sample overlaps by design with the 
THING S and SINGS (Spitzer Infrared Nearby Galaxy 
Survey, Kennicutt et al.l 2003 1 samples and it also has 



behind the disk; the crowding by objects in the fore- 
ground disk and consequently the confusion in the iden- 
tification of the distant galaxies, and, finally, absorp- 
tion of the light from the background galaxies by the in- 
terstellar dust in the foreground disk. Since we are only 
interested in the last one -the dust extinction-, all the 
other factors need to be mitigated and accounted for. 
HST provides the superb resolution to identify many 
distant galaxies, even in the quite crowded fields of 
nearby spiral galaxies. But to fully calibrate for crowd- 
ing and confusion, we developed the "Synthetic Field 
Method" (SFM), in essence a series of artificial galaxy 
counts under the same conditions as the science field 



1) 



galaxies in common with our previous work (Table ((Gonzalez et all (19981 IHolwerda et al.l [2005a I 



To convert the CO (J=2-l) maps to molecular hy- 
drogen surface density maps, we need the conversion 
factor Xco (alternatively denoted as aco)- For the CO 
(J=l-0) line, this is commonly assumed to be 4.4. The 
ratio between the CO(J=1-0) and CO(J=2-l) line is 0.7 
according to the HERACLES observations. To convert 
the CO(J=2-l) map (in K km/s) into molecular sur 
face density , X C o(2-i) = 4.4/0.7 = 6.3M Q /pc 2 (|Leroy 



et al. , 2008 1. The mean values of the CO(J=2-l) surface 



If we identify N galaxies in a field, we need to know 
two quantities to convert this number into a disk opac- 
ity measurement: (1) the number (No) of galaxies we 
would have identified in this field, without any dust 
extinction but under the same crowding and confusion 
conditions, and (2) the dependence (C) of the num- 
ber of galaxies on any increase of dust extinction. The 
disk's opacity in F814W is then expressed as: 

N 



A! = -2.5 C log 



brightness and the molecular hydrogen surface density 
are listed in Table [TJ 



2.3 HST/WFPC2 Images 



No 



(1) 



If the number of identified galaxies behaved exactly 
as photons, the parameter C would be unity. We have 
found it to be close to 1.2 for a typical field, and No 
to depend the surface brightness and granularity of 



the foreground disk (Gonzalez et al. 2003 Holwerda 



The background galax y counts a r e based on HST/WFPC2gg n^gp From Q \ r ^jg^j dist l^^ xy counts 



data, as presented in Holwerda (2005ak and Holwerda 



et al. (2005b). The footprints of the 12 HST/WFPC2 



fields on the integrated H 1 maps of 10 THINGS galax- 
ies are shown in Figure [T] and we only consider these 
areas of the disks. The HST fields are predominantly 



from the Distance Scale Key Project ( Freedman et al. 



2001), and are therefore usually aimed at spiral arms 



in the outer parts of the main disks, in order to fa- 
cilitate the identification of Cepheids. The final driz- 



zled WFPC2 i mages in F8UW and F555W, from |Hoh| 
werda| ( |2005a[ ), can be obtained at |http : / / archive . 



stsci.edu/prepds/sgal/ and the NASA Extragalac- 
tic Databaself] 



3 Disk Opacity from the Number of 
Background Galaxies. 

The central premise of our method to measure disk 
opacity, is that the reduction in the number of distant 
galaxies seen though a foreground spiral galaxy is a rea- 
sonable indication of the transparency of the disk. The 
number of distant galaxies that can be identified is a 
function of several factors: the real number of galaxies 



in the WFPC2 fields, we can obtain both N and C; 
the first from an artificial count of seeded, undimmed, 
distant galaxies, and the second from a series of artifi- 
cial distant galaxy counts with progressive dimming of 
the seeded galaxies. 

Since we cannot know the intrinsic number of dis- 
tant galaxies behind the foreground disk, we treat the 
cosmic variance as a source of uncertainty in N that 
can be estimated from the observed 2-point correla- 
tion function. This typically is of the same order as 
the Poisson error in the opacity measurement Be- 
cause the cosmic variance uncertainty is substantial, 
improvements in the identification of distant galaxies 
barely improve our errors (see also Holwerda 2005b I . 

To test the general SFM results, we have done sev- 
eral checks against other techniques. The results are 
consistent with those obtained from occulting galaxy 



pairs (Holwerda et al. 2005b I , both the results in Domir 



gue 



et al. (2000); White et al. (2000) as well as the later 



opacities found in |Holwerda et al] ( |2007b[ ). The SFM 
results are also consistent with the amount of dust red- 
dening observed for the Cepheids in these fields (the 
majority of which is from the Cepheid Distance Scale 



4 Similar quality products are now also available from the 
archives at STSCI, the High-Level Archive; |www . hla ■ stsci ■ edii| 



J It depends to a degree on the depth of the data. Conserva- 
tively, for this kind of fields, the total error is about 3.5 times 
Poisson l |Gonzalez et al.||2003^ . 
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Key Project Freedman et al. 2001), the dust surface 



densities inferred from the far- infrared SED ( Holwerda 



et al. 2007a discussed below), and the sub- mm fluxes 

below) . Even with 



from KINGFISH observations (j |3.2[ 
HST, the number of identifiable galaxies in a given 
WFPC2 field is relatively small, a fact that results in 
large uncertainties if the field is further segmented for 
its analysis, e.g., sub-divided into arm and inter-arm re- 
gions. To combat the large uncertainties, we combined 
the numbers of background galaxies found in differ- 
ent fields, based on certain characteristics of the fore- 
ground disks, like galactocentric radius, location in the 



arm or inter-arm regions (Holwerda et al. 2005b), sur- 



by a factor of 1.5 or less. Thus, while there is certainly 
a range of dust densities in each field, we are confi- 
dent that the estimate from the above expression is a 
reasonable mean surface density. 

3.2 Herschel-SPIRE Surface Brightness 

Sub-mm data for all our galaxies are available at the 
Herschel Science Archiv^J the majority taken for the 
KINGFISH projeclQ We therefore check the reliabil- 
ity of the SFM as a tracer of the dust surface den- 
sity by directly comparing the surface brightness mea- 
sured by the Spectral and Photometric Imaging RE- 



face brightness ([Holwerda et al.||2005e[ ), or NIR colour ce iver (SPIRE, |Griffin et alj [2010]), onboard Herschel 



(Holwerda et al. 2007c). Because no uniform Hi and 



CO maps were available until now, we compared radial 
H I profiles to our radial opacity profile in | Holwerda] 
et al. ( 2005c I , but this is far from ideal. Now that the 



THINGS and HERACLES maps are available, we can 
compare the average opacity of an HST field to its mean 
H I and H2 surface densities or, alternatively, rank the 
distant galaxies based on the foreground disk's H I col- 
umn density at their position. 

3.1 Dust Surface Densities 

To convert the above opacity of the spiral disk to a dust 
surface density, we assume a smooth surface density 
distribution of the dust (no clumps or fine structure). 
The dust surface density is then: 

= (2) 

K abs 



with K a bs for Johnson I from Draine (|2003), Table 4; 
4.73 x 10 3 



cm g 



. The mean opacity (Asfm) and im- 
plied mean dust surface densities are listed in Table [TJ 
The value for K a b s changes with the types of grain (and 
hence with environment in the disk) and the Draine et 
al. is a value typical for large grains. Variance in K a t, s is 
not unusual depending on the prevailing composition 
of the dust. 

The screen approximation to estimate the surface 
density is common but in fact the dusty ISM is clumped 
and filamentary in nature with a wide range of densi- 
ties and temperatures. Typically, the distant galaxies 



arc seen in gaps between the dusty clouds (Holwerda 



et al. 2007c). The typical value of Ai ~ 1 (Figure 



5j) corresponds to a surface covering factor of 60%, if 
the clouds were completely opaque. In reality, the disk 
opacity is a mix of covering factor and the mean extinc- 
tion of the clouds (on average r c i ud = 0.4 and cloud 
size 60 pc, Holwerda et al. 2007a). We note that our 



mass estimates agree with those from a fit to the Spitzer 



fluxes with the Li & Draine (2001 ) model (to within a 
factor of two Holwerda et al. 2007a Figure [3]). Draine 



et al. ( 2007 ) note that the addition of sub-mm informa- 



to the opacity as measured by the SFM. 

We used the WFPC2 field-of-view as the aperture 
to measure the fluxes at 250, 350, and 500 (im (listed 
in Table [T]), similar to our measurements of the surface 
density in the Hi and CO data (Eland 2.2). These were 



not aperture corrected because of the unique shape 
of the aperture. Figure [2] shows the Herschel surface 
brightnesses versus the SFM opacities for all three wave- 
bands (the 250 and 500 jum. values are the end points of 
the horizontal bars). To convert the flux in a Herschel- 
SPIRE waveband into a dust surface density, one would 
need both a typical dust temperature or a temperature 
distribution and the dust's emissivity. The horizontal 
bars indicate there is a range of mean temperatures in 
these disks. 

There is a linear relation between the Herschel- 
SPIRE surface brightnesses and the SFM opacities. 
The scatter is much less for this relation than between 
the SFM dust surface density values and those inferred 



from far-infrared SED models ( Holwerda et al. , 2007a 
Figure [3]). Hence, we conclude that the SFM opacities 
are a reasonable indicator for mean dust surface den- 
sity. 

As a qualitative check, we compare the dust surface 
densities derived for a subset of the KINGFISH sample 



by Galametz et al. (2012), their Figure Al, to those 



derived above. Typical values mid-disk for the over- 
lap (NGC 3351, NGC 3521 and NGC 3627), where the 
WFPC2 images are located, are ~ O.3M /pc 2 , which 



appear to be typical (i.e., similar to those in Foyle et al. 



2012). These values lie a factor two below the ones im- 
plied by the SFM ( Table|T]), regardless of t he dust emis- 
sivity used in the Galametz et al. (20121 fits but the 



difference is greater for fits where the emissivity is a 
free parameter. We found similar dust surface densities 
from the SED model inlHolwerda et al. (2007a), based 



on the Spitzer fluxes alone (Figure 3 ). Because all these 
models are based on the Draine et al. (2007) model, we 



tion to such a fit may modify the dust mass estimate 



6 |http : //herschel . esac ■ esa. int/Science_Archive ■ shtml| 

7 Key Insights on Nearby Galaxies : a Far-Infrared Survey with 
Herschel, PI. R. Kennicutt, (see als o | Dale et al. , 2012 Gal ametz] 
|et al.|[2012||Kennicutt et al.| [2uTf| [Skibba et al.||2011| 
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Fig. 2: The Herschel/SPIRE 350 /im mean surface 
brightnesses in the WFPC2 field-of-view. Horizontal 
bars mark the 250 (left) and 500 (right) /im fluxes in the 
same field. The width of the horizontal bar is indica- 
tive of the mean temperature of the dust in each disk 
(a wide bar points to higher mean temperature) . Vari- 
ance around the mean surface brightness in each band 
is substantial due to both Poisson noise and structure 
in the galaxy disk. The lowest surface brightness point 
is NGC 3031, the closest galaxy in our sample. This 
field is right on the edge of the ISM disk (Figure [T]) 
and therefore suffers the most from uncertainties due 
to internal structure and aperture correction. 



made a second check using the MAGPHYS SED model 
(magphys). These dust surface density are to a factor 
ten below the SFM or Draine et al. estimates. These 
fits illustrate the importance of the choice of model 
compared to the inclusion of sub-mm data. 

4 H I Column Density and the Number 
of Distant Galaxies 

To improve statistics, one our tactics has been to stack 
the numbers of galaxies in our fields according to a lo- 
cal characteristic (surface brightness, galactocentric ra- 
dius etc.). Here we combine the number of background 
galaxies, both real and artificial, based on the Hi col- 



Fig. 3: The dust surface density inferred by the SED 



model from Draine et al.|(|2007[) based on Spitzer fluxes 



(presented earlier in Holwerda et al. 2007a ) compared 



to those from the SFM. There is at most a factor two 
difference between these, consistent with the lack of 
sub-mm information in these initial fits. Dashed line is 
the line of equality. 



umn density at their respective positions. If there is a 
relation between disk opacity and H I column density 
resolved in the THINGS RO maps, it should show as 
a preference of the real distant galaxies for a specific 
Hi column density, for example for lower values of £hi- 
The artificial galaxies would not prefer any H I column 
density value in particular. The top panel in Figure 
[4] shows the distribution histogram of real (hatched) 
and artificial (solid) galaxies observed, as a function 
of foreground galaxy Hi column density. The bottom 
panel converts the ratio of real and artificial galaxies 
found at an Hi column density into an opacity, using 
equation [l] with C equal to 1.2. The real distant galax- 
ies identified in the HST images do not show a clear 
preference for a certain Hi column density. Their dis- 
tribution is very similar to that of the artificial distant 
background galaxies. As a result, the inferred opacity 
is constant with H 1 column density. In our opinion, 
this lack of a relation can either be: (1) real, point- 
ing to a break-down in the spatial relation between H 1 
and dust on scales of 6" (corresponding to ~ 0.5 kpc 
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HI Column Density (E HJ , 10 20 atoms/em z ) 




10 20 
HI Column Density (2 H] , 10 20 atoms/cm 8 ) 



Fig. 4: Top: histogram of real (hatched) and artificial 
galaxies, N and No respectively, as a function of H I sur- 
face density, Shi- Because all the WFPC2 fields were 
chosen on spiral arms at the edge of the optical disks, 
the range of £hi is limited. Bottom: inferred opacity 
(Aj) as a function of Hi surface density. The dashed 
line is the relation from Bohlin et al. ( 1978 ) for the 



Galactic total (Hi+H 2 ) gas-to-dust ratio. 



in our galaxies); or (2) an artifact of stacking results 
from different fields at various galactocentric radii in 
different foreground galaxies at diverse distances. We 
note, however, that the deviation from the |Bohlin et al.| 
(1978) relation between column density and extinction 
(dashed line in bottom panel) is strongest for the lowest 
Hi column densities, where our statistics are the most 
robust. In our opinion, this points to that one needs to 
compare to the total hydrogen column density, includ- 
ing the molecular component]^] 



5 Average Column Densities and 
Opacity per WFPC2 field 

Our second approach is to compare H I and H 2 column 
densities to disk opacity averaged over each WFPC2 



8 Our fields are usually centered on a spiral arm (to observe 
Cepheids) and this increases the contribution from molecular 
phase. 



field. Table [T] lists the average opacity value for each 
HST/WFPC2 field, and the Hi and H 2 column densi- 
ties averaged over the WFPC2 field-of-view (the foot- 
prints in Figure [j}. The beams of the H I and H 2 obser- 
vations are much smaller than the WFPC2 apertures 
and we expect any aperture correction to the surface 
densities to be small (Table [l|^J Figure [5j left, plots 
the opacity versus H I surface density; there is no clear 
relation between the two, when averaged over the size 
of a WFPC2 field. There are two negative values in 
our present sample [and the entire |Holwerda| ( |2005a| ) 
sample], that are probably due to cosmic variance in 
the number of background galaxies (a background clus- 
ter) . The opacity values and H I surface densities span 
a reasonable range for spiral galaxies. Cui llandre et al. 
(2001) similarly find little relation between reddening 
and number of distant galaxies, on one side, and H I 
column density, on the other. Figure [5j middle panel, 
shows the relation between disk opacity and mean sur- 
face density of H 2 , inferred from the CO observations. 
There are fewer useful points, as there are no CO data 
for three of our WFPC2 fields, and two of the WFPC2 
fields show the aforementioned negative opacity. There 
could be a relation between CO inferred molecular sur- 
face density and opacity. Figure [5j right panel, shows 
the relation between disk opacity and mean surface 
density of total gas (Hi+H 2 ). For those galaxies where 
no CO information as available, we use the Hi mean 
surface density (open diamonds). For comparison, we 
show the canonical Galactic relation from Bohlin et 
al. Opacity appears mostly independent from total gas 
surface density but with the majority of our points lie 
above the Bohlin et al. relation. There is surprisingly 
little of a relation between the gas, total, molecular 
or atomic, and disk opacity. In part this may be due 
in part to the different dust dumpiness in each disk, 
which is observed at a different distance. Alternatively, 
the metallicity and implicitly the average galactocen- 
tric radius of each field is the missing factor in the 
gas-dust relation in these fields. 

One explanation for the lack of a relation in Figure 
[5]is that the measurements were taken at various galac- 
tocentric radii (and hence metallicity) in each disk. Fig- 
ure [6] plots the ratio between the dust surface density 
(to facilitate direct comparison) and the two phases 
of the hydrogen, atomic and molecular in M /pc 2 , as 
a function of radius, scaled to the 25 mag/arcscc 2 B- 
band isophotal radius (R25) from de Vaucouleurs et al. 
(1991). The relation with atomic phase is consistent 
with a constant fraction of Ed /Shi ~ 0.1 with two 
exceptions at R ~ 0.5R 25 ; NGC 3351 and NGC 3627. 
Both of these are small H 1 disks, of which the WFPC2 
field covers a large fraction (see Figure [TJ, both with 

9 We chose not to correct the surface densities because of the 
odd shape of the aperture. Depending on how one treats the 
edges of the aperture, the average surface density varies with 
~10%. 
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Fig. 5: The relation between mean Hi, H2 and total gas (H1+H2) column density (Shi) and average opacity 
(Asfm) f° r each WFPC2 field. Two fields are on average negatively opaque, an effect of cosmic variance in the 
field of galaxies behind them (open circles) and those fields without H2 information denoted by open diamonds. 
There is no clear relation between H 1 surface density and the opacity of the WFPC2 fields, and only a hint of a 
relation between H2 surface density and disk opacity. Most points lie above the the relation found by Bohlin et 
al. (1978), dashed line, right. 



prominent spiral arms. NGC 3627 is a member of the 
Leo triplet and as such may also be a victim of atomic 
gas stripping or a tidally induced strong spiral pat- 
tern. The top panel in [6] shows the ratio between dust 
and molecular surface density, consistent with a con- 
stant fraction of Ed/£h~2 ~ 0.75 with one exception; 
NGC 3198, which is a very flocculant spiral, which a 
much lower Hi surface density. There is little relation 
between the dust-to-atomic or dust-to-molecular ratio 
and radius. Exceptions seem to be either strong spiral 
arm structure, in the case of Hi, or very flocculant, in 
the case of H2 . 

By combining the surface densities of Hi and H 2 
into a single hydrogen surface density (Ehi+ik), we can 
now directly compare the total dust-to-gas surface den- 
sity ratio. In the cases, where no CO observations are 
available (NGC 3031 and NGC 3621), we use the ra- 
tio with H I only. Figure [7] shows the dust-to-total-gas 
ratio as a function of radius. The anomalous ratios of 
NGC 3351, NGC 3627 and NGC 3198 in Figure|6]now 
fall into line. 

If we take the points without CO information (open 
diamond symbols) at face value (assume no molecular 
gas), Figure [7] suggests an exponential decline of dust- 
to-total gas; S d /S H i+H2 = 0.52 x e -4-0-R/H 2 5 , a de- 
cline of the dust-to-total-gas ratio would be consistent 
with the relation with metallicity shown in |Leroy et al.| 
(20111; Sandstrom et al. (20111, and with the trends 



with radius in the recent Spitzer (e.g., Bendo et al 



2010a 



suits (Pohlcn et al. 2010 Smith et al. 20101. 



Munoz-Mateos et al. 2009a I and Herschel re- 



However, if we exclude those points without CO in- 
formation (open diamond symbols) and those with neg- 
ative SFM measurements (open circles), Figure [7] is in 
agreement with a constant gas-to-dust ratio of 0.043 ± 
0.02 (weighted mean). One can reasonably expect a 
much more substantial contribution by the molecular 
component in the inner disk, which would bring the 
three points without CO information into line with this 
constant fraction. This dust-to-total-gas fraction is ap- 
proximately a factor two above the typical value in the 
literature (-0.01-0.03, |Leroy et al.| |2011| | Smith et al. 



2010) or the one from Bohlin et al. (1978). The fact 



that the ratio between dust and total gas surface den- 
sity is nearly constant points to dust in both the diffuse 
H 1 disk as well as in the denser molecular clouds. 



6 Discussion 

When compared to either phase of hydrogen in these 
disks, atomic or molecular, the dust density implied by 
the disk opacity mostly point to a constant ratio. Ex- 
ceptions seem to point to a change in gas phase due to 
the strength of spiral arms in the WFPC2 field-of-view; 
a strong spiral density wave moves gas into the molec- 
ular phase and a flocculant structure into the atomic 
one. A scenario consistent with the density wave origin 
of spiral structure. In our opinion it illustrates the need 
for a constraint on both gas phases for a comparison 
with dust surface density. 

Our dust-to-total-gas ratio of 0.043 (Figure [7]) is 
higher than the values found, for example, in the Lo- 
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Fig. 6: The ratio of dust surface density (£d) to either 
molecular (H2 £h 2 , top panel) or atomic hydrogen (H 1 
Shi) as a function of galactocentric radius, normalized 
to the 25 mag/arcsec 2 B-band isophotal radius (R25 
from de Vaucouleurs et al. 1991 ). Open circles are the 



negative disk opacities from NGC 925 and NGC5194-1. 
The innermost three points from NGC 3031 and NGC 
3621 do not have CO information. 



cal Group spiral galaxies (the Milky Way M31, and 
M33 in the case of Leroy et al.| 2011), or in a single 
Virgo spiral galaxy with Herschel (NGC 4501, Smith 



Fig. 7: The ratio between implied average dust surface 
density (Ed), and the total hydrogen surface density 
(Ehi+h 2 ) as a function of radius. The inner three points 
(open diamonds) do not have CO data. The dashed- 
dotted line is the ratio from Bohlin et al. (1978). 



ity is also a function of cloud cover |Holwerda et al.[ 
2007c). However, our comparison between dust surface 



densities from an SED fit and the number count of dis- 



tant galaxies showed good agreement ( Holwerda et al. 
2007a[ ) (Figure [3]) to within a factor two. We note that 



these SED fits were done without sub-mm information 



et al. 2010 1 or the values found by ?. These studies but Draine et al. (20071 point out that dust masses 



find the values closest to ours in the outskirts of the 
respective galaxy disks. There are several explanations 
for the high dust-to-gas ratio in our measurements: (1) 
we overestimated the dust surface density, (2) a sub- 
stantial aperture correction of the CO and H 1 surface 
densities is needed, (3) for large portions of the disk, 
a different CO-to-H 2 conversion factor (Xco) is appro- 
priate, and (4) a different absorption factor (K a b s ) for a 
disk average is appropriate. First, we are confident that 
our dust surface densities are unbiased and reasonably 
accurate because we checked them agains several other 
observational techniques (Cepheid reddening, occult- 
ing galaxy results, Spitzer FIR SED fits). Our main 
assumption is that the dust is in a screen, which is a 
very rough approximation, especially when the probe 
used is the number of distant objects (i.e, the opac- 



can vary with a factor less than 1.5 if the SED of the 
large grains is done with or without sub-mm informa- 
tion (their Figure 12). The dust surface densities are 
therefore not likely to be overestimated by more than 
a factor two. Our comparison with sub-mm fluxes (Fig- 
ure [2] j |3.2| ) seems to confirm this. The SFM estimate 
of the dust surface density may well be the upper limit 
of dust in these disks. 

Secondly, no aperture correction was applied to the CO 
and Hi surface densities. Because the aperture we use 
to measure the CO and H 1 fluxes is the odd shape of the 
WFPC2 camera's field-of-view (Figure [I]), an aperture 
correction is not straightforward. Yet, we estimate that 
the aperture correction cannot change the reported av- 
erage surface brightnesses sufficiently, as the resolution 
of the observations is substantially smaller than the 
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WFPC2 aperture (Table [TJ. 

Thirdly, when averaged over a large portion of the disk, 
which spans a range in density environments, the CO 
conversion factor (Xco) may well underestimate the 
total molecular hydrogen surface density, since some 
molecular clouds the observed CO may be from the 
"skin" of the CMC and there is not straightforward 



conversion from CO to H2 volume ( 


Feldmann et al. 


2011a Glover & Mac Low 2011 


Mac Low & Glover 


2012; Madden et al.[ 2011||Planck Collaboration et al. 


20TTlphetty et al. 2011| |Wolfire et al. 


2010 ?). 



A fourth option is that the dust absorption factor (n a t, s ) 
is different when averaged over different environments 



and therefore dust grain properties (e.g., Narayanan 



et al. 2011 2012), although this is likely a secondary 
effect. 

If all our dust surface densities would are all over- 
estimated by a factor ~2, or the aperture correction in- 
creased the gas surface densities substantially, one may 
not need to change the Xco factor to bring our dust- 
to-gas ratio in line with recent results from Herschel. 
We suspect, however, that the explanation includes a 
different Xco , when averaged of a large section of the 
spiral disk and at different galactocentric radii, extend- 
ing the range in Xco values found in Local Group spiral 




12 + log(0/H) 



galaxies by Leroy et al. (2011) 



6.1 Comparison to Metallicities 

The present consensus is that the dust-to-total-gas ra- 
tio depends linearly on the metallicity (see for instance 



Leroy et al. 2011). Fortunately, uniformly determined 

T7T 



metallicity gradients for the SLNGS^J and h ence THING S 
and HERACLES, galaxies are presented in |Moustakas| 
et al.l (120101). Only NGC 3627 does not have metallicity 



Fig. 8: The logarithm of the ratio between the total 
hydrogen surface density (Shi+h 2 ) and the implied av- 
erage dust surface density (Ed) as a function of the 
metallicity (12 + log(0/H)), estimated from Figure 7 in 
Moustakas et al. ( |2010 ). The circles and diamonds are 
for the calibration from Kobulnicky & Kewley (2004) 



and Pilyugin & Thuan ( 2005 ) respectively. Open sym- 



information. Starting from their linear relation for the 
radial dependence of metallicity in each galaxy (their 
Figure 7) , we can obtain an estimate of the metallicity 
for each of our WFPC2 fields. They present two dif- 
ferent estimates of metallicity (log(0/H)), with cither 
the theoretical calibration from |Kobulnicky fe Kewley 



bols are those points without CO information. There is 
no metallicity estimate for NGC 3627. We use the gas- 
to-dust ratio here in order to compare to the relation 



from Leroy et al. (2011) 



(20041 or the empirical one from Pilyugin & Thuan 
( 2005 1 (see Table [T|). p"oustakas et al.| ( |2010[ ) note that, 



two calibrations. Our points lie lower than the linear 



relation from Leroy et al. (20111 for the gas-to-dust 



until the calibration issues are resolved, one should ei- 
ther average the metallicity estimates based on either 
calibration or use both separately. We will use both 
calibrations separately for comparison and the total- 
gas-to-dust ratio to facilitate a direct comparison to 



ratio with metallicity, not unexpectedly as we already 
established that our dust-to-gas values are higher than 
those previously reported. 

However, using the calibration from Pilyugin & Thuan 



(2005), and discarding those points that lack CO infor- 



mation, there is a reasonable agreement with the rela- 



Figure 6 in |Leroy et al.| ( |2011[ ). We note that since tion from|Leroy et al.| ([20TT I 



our WFPC2 fields were placed with crowding issues in 
mind, our coverage of galactocentric radii (and hence 
metallicities) is not very large. 

Figure [8] shows the logarithm of the total-gas-to- 
dust ratio as a function of metallicity, using cither of the 



10 Spitzer Infrared Nearby Galaxy Survey (Kennicutt et al. 
[2003] ). 



7 Conclusions 

To conclude our "Opacity of spiral disks" project, we 
have compared the opacity of spiral galaxies, and the 
hence dust surface density to the surface densities of 
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hydrogen, both atomic and molecular, the original goal 
of our project. We conclude from this comparison: 

1. The disk opacity scales with the Herschel-SPIRE 
250 /im surface brightness (Figure [2j), confirming 
our assertion that opacity scales with dust surface 
density to first order. 

2. There is little relation between the Hi column den- 
sity and where a distant galaxy was identified in 
these fields (Figure [I]). 

3. Averaged over a WFPC2 field, there is only a weak 
link between disk opacity (or dust surface density) 
and gas surface density, either atomic, molecular or 
total (Figure [5]), pointing to third factor; radius or 
metallicity. 

4. The dust-to-H 1 or dust-to-H2 relations with galac- 
tocentric radius are both relatively constant (Fig- 
ure [6]) , but the exceptions point to the role of spiral 
structure in the dominant gas phase of the ISM. 

5. The dust-to- total-gas ratio is close to constant for all 
our fields S H i+h 2 = 0.043 ± 0.024 (Figure 0. This 
higher value can, in our opinion, be attributed to a 
different conversion to dust surface density or the 
CO-to-H 2 conversion factor (Xco) f° r such large 
sections of disks. 

Compared to the relation between total-gas-to-dust 



(> 



and metallicity from Leroy et al. (2011 1, our results 



are reasonably consistent, provided one uses the|PT 



lyugin & Thuan ( 2005 ) calibration of the metallici- 
ties of |Moustakas et al.| ( |2010[ ) (Figure [§). 



Future use of the number of distant galaxies identi- 
fied through a foreground spiral disk as a probe of dust 



is critically limited by cosmic variance ( Gonzalez et al 



2003 Holwerda et al. 2005d) but its optimal applica- 



tion will be on a single large HST mosaic of a nearby 
face-on spiral (e.g., M81 or M101), which will most 
likely the last contribution of this unique approach to 
the issue of the dust content of spiral disks. 
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A MAGPHYS SED Model 



As an alternative check of the inferred dust masses, 
we ran the Multi-wavelength Analysis of Galaxy Phys- 
ical Properties (magphys) package on the Spitzer and 
Herschel/SPIRE surface brightnesses. This is a self- 
contained, user-friendly model package to interpret ob- 
served spectral energy distributions of galaxies in terms 
of galaxy-wide physical parameters pertaining to the 
stars and the interstellar medium, following the ap- 
proach described in da Cunha et al. ( |2008[ ). Figure 
Al| summarizes the result: dust surface density derived 
from the magphys fit compared to those inferred from 



the number of distant galaxies. In Holwerda et al. ( 2007a I, 



we found that the Draine et al. (2007) model inferred 



similar dust optical depths for these disks as the SFM 
as well as similar (to within a factor two) dust masses. 
The discrepancy with magphys illustrates, in our view, 
the importance of modeling sections of spiral disks in 
resolved observations with more physical models that 
include a range of stellar heating parameters (e.g. the 
models by Draine et al. 2007 Galliano et al. 2011). 




0.05 0.1 0.15 

E a (M s / pc 2 ) - MAGPHYS 



Fig. Al: The dust surface densities from the MAGPHYS 
fit and inferred from the number of identified back- 
ground galaxies (SFM) for each WFPC2 aperture. The 
dashed line denotes a factor ten ratio, magphys SED 
models do not take internal structure and differential 
stellar heating into account. 
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SED of each WFPC2 field with the magphys fit. 
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